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We re-examine results of the Liouville theory and provide arguments that a negative bare cosmo- 
logical constant is essential to define two-dimensional quantum gravity. From this we are naturally 
led to a regularization of quantum gravity within the Regge approach such that it is described by 
small fluctuations around equilateral triangles, whose average link length approaches zero in the 
continuum limit. We investigate a model based on this idea numerically and present evidence for 
the desired long-range correlations. Interestingly, the approach might generalize to higher dimen- 
sions. The picture of an inflated balloon, which is often used to demonstrate the properties of an 
expanding classical universe, seems to be valuable to understand quantum gravity as well. 

PACS: 04.60-m. 



I. INTRODUCTION 

In recent years, substantial progress has been achieved 
in understanding two-dimensional quantum gravity 
based on the Liouville theory jj]], see Q for a review. 
Unfortunately it was not possible to extend these results 
to higher dimensions. In the light of these difficulties we 
re-examine the results of Polyakov, starting with the path 
integral of two-dimensional gravity for a fixed topology 
and withf] D additional matter fields X 1 : 



Z D {lM)) = J VgVX 1 



exp 



J d 2 x^ - I d 2 x^g ab d a X l d b X" 
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Following the procedure of Polyakov, we 'fix a gauge' us- 
ing coordinates x so that g a b — e^h a b, where h is an 
external background metric. Using an appropriate regu- 
larization, this leads to the well-known path integral 



Z D (p a ) = / V<j>{x) x 
exp 



(2) 



for the Liouville field 4>(x), with R being the curvature 
of the metric h, which is chosen so that R is a constant 
depending on the topology JIJ]. The conformal anomaly 
is present for D ^ 26 where a scalar field appears in the 
quantum theory although there are no physical degrees 
of freedom in the classical theory. (There are no field 
equations for gravity in two dimensions.) The parameter 



'Following established notation [g], in the present context 
D denotes the number of matter fields or dimension of the 
embedding space. 
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where e is the cut-off parameter of the regularization 
procedure. Usually the path integral (||) is the start- 
ing point for further examination. We emphasize only 
that for R < the interaction term R<fi + /j, s e^ has a true 
minimum at 4>$ = ln(— R/ fJ. s ) and using 4>{x) — 4>o + r)(x) 
one can approximate the vacuum functional by 



exp 
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with c = -\R etc. g. 

We highlight the fact that for D < 2 relation (||) indi- 
cates that the bare cosmological constant [1q is negative 
and infinite in the limit e — *• 0. While the path inte- 
gral (Q) seems to be highly divergent for /io — > — oo, the 
Polyakov procedure transforms it so that the scalar Li- 
ouville field results in two dimensions. In this paper we 
want to understand how this is possible, using the frame- 
work of the Regge calculus || . Subsequently, we touch 
briefly on the prospects of generalizing our regularization 
procedure to higher dimensions. For our investigations 
of quantum gravity in the Regge regularization we use 
a simple measure, thus defining a model. Although it is 
not clear whether this model can be fully equivalent to 
the Liouville theory in two dimensions, it shares some of 
its remarkable properties and exhibits a phase transition 
which appears to define a continuum theory. However, 
it might be that two dimensional gravity is a poor test- 
ing ground for the general case |@|. In higher dimensions 
the Einstein Hilbert action should dominate the physics 
of quantum gravity and since, without proper regulariza- 
tion, it is unbounded in the Euclidean sector, its behavior 



1 



needs to be determined before one can discuss the sub- 
tleties of the correct measure. 

In the next section we define the model and obtain ana- 
lytical results including properties of its phase transition. 
Numerical simulations are subsequently performed in sec- 
tion III . They yield long-range correlations as needed for 
the continuum limit. Conclusions follow in section IV. 



II. THE MODEL 

For D — we can write the path integral ([!]) as 

Z{ l x )= [ Dge-^ (5) 



where A [g] is the total area of the geometry. The inte- 
gral (^|) is defined only after proper regularization as the 
following example Q illustrates: Consider a smooth ge- 
ometry of finite total area. Then, assume that a "spike" 
of length L and circumference u grows out of this area. 
The spike can have arbitrarily large length L, but, if u 
is sufficiently small, it would not be suppressed by the 
action of the path integral due to its infinitesimal area 
~ uL. Indeed, the path integral would be dominated by 
ill-defined geometries, unless we introduce a regulariza- 
tion procedure, which provides for a cutoff to prohibit 
L from becoming arbitrarily large. A proper regulariza- 
tion procedure introduces counter-terms, which prevent 
degenerate geometries and lead to a well-defined path in- 
tegral Q , as we will discuss in the following for a model 
of pure gravity. 

We consider a triangulation of fixed topology (e.g a 2- 
torus), with N2 triangles and N\ links. We assume that 
the triangulation is sufficiently regular and, following the 
Regge approach, consider the link lengths x as the vari- 
ables of this model. We use a cutoff a on the link lengths 
and the simple local measure Yli dxi to obtain the path 
integral 
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where the total area A[x] of the geometry is the sum of 
all triangle areas, 



and the integration range restricts the link lengths to 
< x < a. The above integral is obviously well-defined 
for all /io as long as a is finite. 

It is not clear whether the partition function (|^) can 
be in the same universality class as (||), i.e. that dif- 
feomorphism invariance will be properly implemented by 
our simple measure in a continuum limit. It is known, 
however, that the Regge approach approximates Einstein 
gravity in the classical limit || . We reach the continuum 



limit formally by increasing the number of links N\ and 
decreasing their lengths, a — > 0. 

To proceed further, we re-write the path integral (||) 

as 



Z(jio,a,n) 



/ n dxi exp _ a*o^n 



1 



(xi/ a y 



(7) 

i.e. replacing the cutoff with a counter-term, which is 
equivalent in the limit n — > 00. 

Let us qualitatively examine how the Regge lattice 
might behave for different values of /io. A positive cos- 
mological constant, /io > 0, tends to suppress large ar- 
eas, while the counter-term restricts the link lengths. In 
the limit fio — > +00 the expectation value of the total 
area will tend to 0, but it is not clear whether the link 
lengths will approach as well. It is possible that their 
expectation value remains finite and this would result in 
a 'crumpled' lattice with collapsed triangles. The numer- 
ical simulations presented in the following section suggest 
that this indeed happens and it is not clear if a reason- 
able continuum limit can be found for /io > 0. Previous 
numerical investigations Q have been mostly confined to 
this range of a positive cosmological constant. 

On the other hand, if /io < 0, the cosmological term 
tends to inflate the area of the lattice, while the counter- 
term limits the link lengths. Therefore, equilateral trian- 
gles are preferred and for /io — ► — 00 we have to expect 
a lattice which consists of (almost) equilateral triangles, 
with small fluctuations around the maximum link length 
a; other configurations cannot significantly contribute to 
the path integral. The situation is thus similar to an 
inflated balloon, where the pressure of the air is the ana- 
logue to the cosmological term, while the role of the rub- 
ber molecules is played by the Regge links. If we inflate 
the balloon enough, the rubber becomes locally flat and 
elastic. While it is immediately clear that the Regge lat- 
tice is locally flat for regular triangulations in two dimen- 
sions fIo| , if /io is sufficiently negative, we must examine 
in the following whether the lattice fluctuations are able 
to reproduce a scalar field as in the Liouville theory. To 
this purpose we consider first the path integral for /io < 
and n > 2. The special case of n = 2 is then discussed 
separately. 



A. n > 2 

If n is large but finite and /io is sufficiently negative, 
only lattice configurations with link lengths x ~ a con- 
tribute significantly to the path integral. Configurations 
with x » a are suppressed by the term —(x/a) n , while 
configurations with x « a cannot contribute signifi- 
cantly due to the area term — /to A > 0. We thus replace 
xi with a(l+£i) where the variable determines the fluc- 
tuation of link / around the length a. The path integral 
becomes 
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Z(po,a,n) = a Nl / JJc^i x 
J i 

exp (^ a 2 £ A t [l + £] - E(l + 



(8) 



and only configurations £ ~ contribute significantly. 
Notice that the cut-off a appears as a simple, unphysical 
multiplicator of the bare cosmological constant, which 
allows us to remove it from the path integral by rescaling 
Mo- 

The area of a single triangle with the links a, 6, c is of 
the form 

A t [l + £] = A a + Ai(f a + & + &) - A 2 (C + $ + 

+ (6 - Cc) 2 ) + o(e) (9) 

with positive constants 
1 
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The cosmological constant term, i.e. the total area of the 
geometry, can therefore be expanded as 

E A t [l + e] = A N 2 + 2Ai E& - 2A 2 E £ 2 

-s 2 ^(^-e«) 2 + o(e 3 ) (ii) 



where the sum over [nm] indicates a sum over neighbor- 
ing links. Correspondingly we expand the regulating link 
term as 



n(n — 1) 



E^ 2 + °« 3 ) 



(12) 



and insert both expansions in the path integral. We as- 
sume equilibrium of the link lengths around a and for 
consistency we have to set u a 2 2Ai = —n, so that the 
first order term in the action vanishes. This determines 
the bare cosmological constant as 
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so that it is negative and tends to — oo for a — > 0. (Notice 
the similarity with relation ([|).) 
The path integral becomes 



Z(n) = const / d£i exp 



n(n-l) A 2 



B 



E^ 2 

i 

(14) 



where we now understand n as a coupling parameter 
which is not necessarily an integer. While the first term 
guarantees that |£| <C 1 for large enough n, the second 
term can obviously be interpreted as the discretization of 
the kinetic term of a scalar field £, similarly to the field 
i] of the path integral (H|). 

To proceed further, we discretize the variables using 
eo~i instead, with o~i being an integer ±1, ±3, ±5, so 
that the difference between two neigbouring values of £ 
equals 2e. The path integral is then replaced by a sum 
over different configurations o~\ 



z(/3, e ) = E ex pf-^ 2 E CT2 
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where we use (3 = ne 2 and assume that n is large enough 
that the coupling of the first term in (|l4) is essentially 
n 2 /2. 

The limit e — * leads back to the path integral (|l4| ) 
and, if we adjust n such that (3 remains finite, the (3 2 e~ 2 
term of ([l5|) will suppress all higher values of a;. It is 
then sufficient to perform a summation of all configura- 
tions a i = +1,-1 and the system becomes equivalent 
to an Ising model. Therefore, in this limit the model 
exhibits a second order phase transition at a certain crit- 
ical coupling /3 C , which leads us to the equivalence with 
an interacting scalar field with 5*-matrix —1 H]. Al- 
though this field is in a different universality class than 
the Liouville theory, the associated long-range correla- 
tions demonstrate that our 'quantum balloon' is indeed 
elastic as we had hoped for. In previous investigations it 
has been proposed to approximate quantum gravity by 
Ising models and numerical simulations of such models 
have been performed [fl2|| . Our calculation demonstrates 
that these Ising models follow naturally in a certain limit 
(po — > — oo, n — > oo) of the Regge approach. 

The occurance of different universality classes is cer- 
tainly related to the sensitive dependence of two dimen- 
sional quantum gravity on the definition of the measure. 
The above calculation is supposed to depend on the de- 
tails of the limit (fi — > — oo,n — > oo) and, furthermore, 
we also do not know the universality classes for general 
n. In principle, this could be determined by calculating 
the critical exponents numerically. So far, we did not in- 
vestigate this issue in detail. It is expected to be of minor 
importance in higher dimensions due to the presence and 
dominance of the Einstein-Hubert term there. 



B. n = 2 

Having discussed the case n > 2 and a particular limit 
n — > oo, we examine the special case n = 2 in the fol- 
lowing. The reason for our interest in this case is that 
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the critical point is found for a finite value of fio, which 
allows for easier numerical investigations. 
We use the dimensionless variables 



qi = (xi/a) 2 



(16) 



with the measure Yii dqi instead of Y\i dx\ for reasons of 
analytical simplicity. Using the rescaled coupling param- 
eter 



the path integral (0) becomes 



Z(jj) 



\ [ dqi exp -fiA[q] - ^ q t 



(17) 



(18) 



This path integral is well defined as long as \i > ji c and 
the value of the critical coupling fj, c is found from the con- 
figuration with all link lengths equal, qi = q. The action 
reduces to (—[i(y3/2)N2 — N%)q in this case and the ex- 
ponent remains finite only if /i > \i c = —(2/y3)(Ni/N2)- 
In two dimensions we have N\ = 3N2/2 and 



-V3 



(19) 



is independent of the lattice size. Although the lattice 
is expected to fluctuate heavily at the critical point, we 
have nevertheless determined \i c by assuming that config- 
urations around equilateral triangles dominate the path 
integral. Unfortunately, we cannot prove this assump- 
tion, but it is confirmed by our numerical simulations 
(see the next section) and it is self-consistent as we show 
in the following. 

We substitute the N\ variables \q\, gjvj by the Ni 
variables [<7, £1, ...,£_/Vi-i] defined by 



Nx-l 

and q Nl =q(l- £;) , 
1=1 



(20) 



(21) 



so that J2i Ql — N\q. The total area A[q] can then be 
written as 



A[q] = (VS/2)N 2 qr[^ 



(22) 



where the function r depends on the variables only 
and varies between and 1. It equals 1 if all triangles 
arc equilateral (maximum area) and it is if all triangles 
collapse (minimum area). The path integral becomes 



J Y[dtiexp(-^(iN 2 qT[Z}-Niq 



(23) 



because the functional determinant is a constant, when it 
is evaluated under the assumption of small fluctuations 
around equilateral triangles, which allows us to ignore 
the triangle constraints. Integrating out q, one arrives at 
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We are interested to determine the expectation value of 

(25) 



2 N x 

<r> for M ^ Mc = -__ 



and obviously the behavior of the function r[£] near r = 1 
is essential in this case. We use the fact that r has a 
maximum at £ = and expand it as 

r[£] = 1 - BN^ 1 ~ E B ^ + °^ 3 ) 

1 M 

with a constant B > and a constant matrix B^j. In 
a second step we introduce spherical coordinates £j = 
rsin^i sin02--- and thus arrive at 



Z(fi) — const J Y\_d<fidr 



r N '- 2 x 



(26) 



J{6) 



AT, 



if we neglect terms of order 0(£ 3 ). The function J(4>) 
denotes the Jacobian for the spherical coordinates and 
fij denotes a matrix of functions which depend on the 
angles 4> only such that = r 2 fij((j)). At fi = fx c the 
path integral becomes 



Z(/i) = const J J^J i 



J(<j>) 



(B + EwBijfiM 

r Nx-2 



dr ■ 



r 2N t 



(27) 



which we can write as const J dr r -( N i+ 2 ) [ n order to 
calculate the expectation value < r >. Obviously the 
path integral has a strong singularity at r = and the 
expectation value of r is zero at fi c . This, of course, 
implies that < r >= 1. 

In order to determine < r > for /i = \i c + S, the first 
momentum J drf(r) rj J drf(r) needs to be calculated 
for the function 

f{r)=r N ^ 1 /(S + B'r 2 ) N ' , 

where B' equals B + J2[ij] fij ["!>]• The function / has a 
maximum at r max = y/(N~i - VjEJiN^B 1 + W) and falls 
off sharply away from it. We conclude that < r >^ 
const 5 2 where the constant depends on the integration 
of B' over the spherical angles cf>. Therefore, 
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1 — const (/i — /i c )" 



(28) 



for lattice configurations fluctuating around equilateral 
triangles, providing for the self-consistency of our initial 
assumption. 

An eventual, non-trivial continuum limit fi — > fi c is 
quite subtle. The numerical investigations of section |T| 
suggest the following picture: At \i c the path integral is 
dominated by small fluctuations around equilateral trian- 
gles. The link lengths show correlations which for /i — > u c 
are consistent with a diverging correlation length (mass 
to zero). Considering the ratio r = mm(ir;)/ max(xi), one 
hnds that the expectation value < r > shows a strong fi- 
nite size effect, in contrast to < t > which is essentially 
independent of N\. The continuum limit is then found 
for (Ni — > oo, n — ► (a c ), so that < r >^ 1 and < r >^ 0, 
which indicates that the triangles are (almost) equilateral 
yet over large distances the geometry of a typical lattice 
is non-trivial. 



C. Higher dimensions 

Let us briefly comment on the situation in higher di- 
mensions. We focus on the physically interesting case of 
four dimensions and, because of the importance of nu- 
merical calculations, on n — 2. We consider the path 
integral 



Z(n) 



J Dq exp ( -nV[q] 



(29) 



with total 4- volume V[q] and a regulating counter-term. 
We use again a simple, local measure for simplicity and 
we ignore the Einstcin-Rcggc action as well as other pos- 
sible terms in the action. We introduce variables p and £ 
so that pNi — J^i if an d q% — P 1 (1 + &) and we use the 
function v[£] as we used r for two dimensions, so that 
Viol ~ = P^K^oA/4 where the constant V a is 

chosen so that the maximum of v equals f. The path 
integral becomes 



Z(ji) 



dpp 



(N 1 -2)/2 



J ex p (-^v p^] - n 1P ) 



(30) 



Proceeding along the same lines as for the two- 
dimensional case, one concludes that < v >^ 1 as u 
approaches the critical value u c . At the critical point lat- 
tice configurations with equilateral 4-simplices contribute 
significantly only, indicating long-range correlations and 
a non-trivial phase transition as in two dimensions. 

Unfortunately, a regular triangulation with equilateral 
simplices in four dimensions does not provide for a locally 
fiat geometry. Also, the ratio N4/N1 is not constant and 
some properties of our regularization depend on the de- 
tails of the triangulation. Clearly, a mechanism is needed 



which would favor locally flat geometries, independent of 
the triangulation and we hope that the Einstein- Hilbert 
action will achieve this. Indeed we can show that our 
regularization is sensitive to the Einstein-Hubert action 
near the critical point, by adding 



(31) 



to the action in equation (29). Using the gravitational 
coupling (3 we write 



(3S EH {q) = -/3p 1/2 s(0 



(32) 



where s(£) is bounded because the deficit angles —nit < 
5t < 2tt are bounded (n depends on the details of the 
triangulation) for a given triangulation. 
One can expand the action for small (3 as 



exp(-/V/ 2 S (0) - 1 - P P lrz s{i) + O(p) 



1/2. 



and evaluate the path integral equation ( |30[ ) with this 
correction. Performing the integration of the variable p 
one finds that the integrand in the remaining path inte- 
gral is modified by a factor 

[l-/3C(JVi) *&)]/[!- (/V/OK6)] 

where C(N\) is a constant, which depends slightly on Ni. 
Near the critical point we can set £ = to approximate 
(3 C as the gravitational coupling which leads to 

[ 1 - & C(Nt) 8(0) } I [ 1 - (/x/Mc) i/(0) ] = . 

Of course, this does not imply that Z = at /3 C , rather 
we have to conclude that our approximation for small (3 
breaks down. It is now important to see that /3 C — > 
for fi — > /x c , which means that our model becomes very 
sensitive to the Einstein-Hilbert term near the critical 
point. Numerical simulations could clarify the relevance 
of this effect and in general the phase structure for the 
two coupling parameters \i and (3. This is left as a project 
for future investigations. 



III. NUMERICAL RESULTS 

To support the assumptions as well as the results of 
our analytical calculations in the previous section, we 
perform numerical simulations of our two dimensional 
model. We use the path integral (0) for n = 2 on a 2- 
torus. Our consideration which leads to u c = — v3 of 
equation ( |l^ ) is not affected by the change of the mea- 
sure from the path integral Jig ) to (0). We have gen- 
erated data on 8 x 8, 16 x 16 and 32 x 32 lattices at 
fj, = 5.0, 2.0, 1.00.0 - 1.0 - 1.5 - 1.7 and -1.72. Each 
data point relies on at least 1M sweeps and measurements 
are taken every 10 sweeps to reduce auto-correlation of 
the data. 
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FIG. 1. Together with 2{A) / (N 2 V3{xi)^) the ratio 
(r) = 2{A/xf }/(A r 2V / 3), defined by equation (p2j), is depicted 
as function of the cosmological constant /i. Within the accu- 
racy of the plots data from 8 x 8, 16 x 16 and 32 X 32 lattices 
fall on top of one another. 
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FIG. 2. The expectation value of the ratio 
r = min; (xi)/ max; (xi), where minimum and maximum are 
defined with respect to the configuration at hand, is plotted 
versus fj, for lattice sizes 8 x 8, 16 x 16 and 32 x 32. 



In figure |l| the expectation value 

(r) = 2(A/xf)/{N 2 V3) (lower curve) 

is shown together with 

2(A)/(N 2 V3(xi} 2 ) (upper curve) 

for various values of the cosmological constant. There are 
almost no finite size effects for these quantities. Within 
the accuracy of the plot our data from 8 x 8, 16 x 16 and 
32 x 32 lattices fall on top of one another. The curve 
for (t) decreases for positive values and tends towards 
zero for increasing bare cosmological constant (fx — > oo), 
indicating a crumpled lattice consisting of collapsed tri- 
angles. The ratio is finite for negative values of the cos- 
mological constant and approaches for fi — > fx c = — Vo its 
maximum possible value 1, indicating the dominance of 
equilateral triangles in a locally flat lattice as predicted 
by equations (|l^) and fl2q). Link length fluctuation are 
demonstrated by the fact that the curves for (r) and 
2 ( A)/ (N 2 V3(xi) 2 ) are quite distinct. 

On a single configuration the range of link fluctua- 
tions is between x^ m - ln = mhii(xi) and xi >max = max;(x;), 
which define the minimum and maximum link length on 
the configuration. In figure ^| we plot the average value 

— (^^,min/^'Z,max) 

for our three lattice sizes. As we assumed, the expec- 
tation value of this ratio decreases with increasing lat- 
tice size and the results show that even a modest lattice 
size (32x32) is sufficient for substantial variations of the 
overall geometry near the critical point ((r) < 0.1 for 
(r) > 0.9 ). 



We are now confident that the assumptions of the pre- 
vious sections are essentially correct, leading to the pic- 
ture of a lattice with highly correlated links near the crit- 
ical point. In the following we check this conclusion, in- 
vestigating the expected long-range correlations directly. 
We employ that we know the exact value (|ig| ) of fi c to 
control the approach to the critical point. To improve 
the numerical stability we add a n = 3 term fxf and 
investigate the partition function 



ZQ*,f) = 



I T\dxi x 
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exp 



A, 
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(33) 



We take fi = \i c and are interested to establish critical 
behavior for / — ► 0. As long as / > 0, the integral ([53]) is 
well-defined for all /i. However, for / = the lattice has a 



finite area for fx > fi c only, as discussed in subsection p B . 

We compute the link-link correlation functions defined 
by 



cm = 



((x xt) - (x )(x t )) 

((xqXq) - (x Q )(x )) 



(34) 



where, using the regular triangulation [[10|, two face to 
face links pointing in the x direction and separated by 
t lattice spacings in the t direction are correlated. Note 
that t is not necessarily proportional to the geometrical 
distance between the links, especially for large t, since 
we know that the geometry fluctuates according to our 
measurement of (r). 

Figure || compares our 32 x 32 and 64 x 64 results and 
shows indeed a strong correlation of the link lengths for 
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FIG. 3. The 2-point functions @ on 32 x 32 (t < 16) and 
64 x 64 (t < 32) lattices are plotted for / = 1CT 7 , 1CT 6 and 
1(T 5 . 



FIG. 4. The 2-point functions @ on a 32 x 32 lattice are 
plotted for our higher / values, / = 10~ 4 , 10~ 3 , 10" 2 and 
10" 1 . 



/ — > 0, even for large t. However, irregular finite size 
effects are present as well (e.g. for / = 10~ 6 ) and a 
full-scale FSS analysis, which is beyond the scope of this 
paper, seems necessary to clarify the critical behavior of 
our system. 

For the 32 x 32 lattice we performed also calculation 
at higher values of /. The long distance correlations dis- 
appear then quickly, as is we show in figure Hfnote the 
change in ordinate scale between figure |] and|]) . 

In order to determine the mass of the 'particle' asso- 
ciated with the lattice fluctuations, we extracted mass 
values using so called cosh fits 



Ci(t) = A (e- ,nt + e-™^-^ 



(35) 



from the correlation functions. Figure |^ depicts our re- 
sults from the 32 x 32 and 64 x 64 lattices in the limit 
/ — *■ 0, including also / = 10~ 4 from our larger / values 
on the 32 x 32 lattice. The mass values agree down to 
/ = 1CP 6 ; for / = 1CP 7 we find then disagreement, which 
is expected when the small lattice can now longer accom- 
modate the correlation length. Again, we emphasize that 
our data are currently not sufficient for a complete FSS 
analysis of the model. 

After we have demonstrated numerical results in two 
dimensions, we are ready to have a peep at higher di- 
mensions. Figure (|^) shows very preliminary results 
of numerical calculations in four dimensions. We per- 
formed these calculations without gravitational coupling 



(i.e. P = 0) near \x 



and include a term e qf in 



the partition function (E9h for improving convergence. 

The figure shows the results with error bars from two 
independent simulation runs ona4x4x4xl0 lattice, 
measuring the correlations for the long direction. Each 
simulation performed 10k sweeps through the lattice and 



sampled every 10th configuration. Depicted is the corre- 
lation function 



CM) 



((ffogf) - (go)(gt)) 
((gogo) - (go) (go)) 



(36) 



for the links forming the edges of the hypercubes in the 
lattice. The correlations are of course measured in the 
long direction. Although the figure suffers from large 
statistical errors, it exhibits the desired long-range cor- 
relations. 



IV. CONCLUSIONS 

We defined a class of models of fluctuating geome- 
tries, based on the Regge approach, to illustrate basic 
properties of the Polyakov path integral in two dimen- 
sions. Although we used a simple, local measure we find 
evidence for the existence of a non-trivial phase transi- 
tion with long-range correlations and the emergence of a 
scalar field. We demonstrated the importance of a nega- 
tive bare cosmological constant, which may open the door 
to investigations of related systems in higher dimensions. 
Our calculations indicate that one could expect a non- 
trivial phase transition with long-range correlations in 
four dimensions as well and numerical simulations with 
negative cosmological constant in higher dimensions are 
the obvious next step to investigate these models further. 
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FIG. 5. The mass fits @ from 32 x 32 and 64 x 64 lattices 
in the limit / — > 0. 
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FIG. 6. The 2-point functions @ is shown on a 4 3 10 lat- 
tice for (5 = 0, fj, = —26.5 and e = 0.01, where the correlations 
are measured along the long direction. 



The simulations were done on workstations of the FSU 
HEP group. 
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